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ABSTRA CT 

A  scale-invariant  theory  of  Raman  scattering  of  light  by  fractal  clusters  is  de¬ 
veloped.  The  enhancement  factor  GKS  of  Raman  scattering  is  shown  to  scale  in  terms 
of  a  properly  chosen  spectral  variable  X.  The  critical  indices  of  the  enhancement  factor 
are  found  to  be  determined  by  the  optical  spectral  dimension  of  the  fractal.  Numerical 
modeling  is  carried  out  and  shown  to  support  the  analytical  results  obtained.  The  the¬ 
ory,  which  does  not  contain  any  adjustable  parameters,  agrees  well  with  experimental 
data  on  sutface-enhanced  Reiman  scattering  over  a  wide  spectral  range. 


u'A 


PACS:  78.30.Ly,  61.42+h,  78.65.Ez,  73.20.Mf 


1.  INTRODUCTION 


Surface-enhanced  Raman  scattering  (SERS)  is  one  of  the  most  intriguing  optical 
effects  discovered1,2  in  the  past  15  years.  Its  manifestation  is  a  very  great  increase,  often 
as  large  as  a  million-fold,  in  the  Raman  scattering  intensity  of  molecules  absorbed  on 
the  surfaces  of  special  metallic  structures,  as  compared  to  the  intensity  of  an  equivalent 
number  of  molecules  in  solution  or  the  gas  phase.  The  special  metallic  structures  that 
display  SERS  normally  involve  small  metal  features  such  as  those  found  on  micro¬ 
rough  metal  surfaces  and  in  aggregated  colloidal  particles,  whose  dimensions  are  much 
smaller  than  wavelength  A.  The  major  contribution  to  the  enhancement  is  understood 
to  originate  from  the  large  local  electromagnetic  fields  that  arise  from  resonant  optical 
excitation  of  surface  plasmons.  Because  these  plasmons  are  both  intense  (i.e.  possessing 
great  dipole  moments  and  low  decay  rates)  and  easily  excitable  with  commonly  available 
laser  sources  for  the  alkali  and  coinage  metals,  the  SERS  effect  is  largely  restricted  to 
these  metals. 

Theory  shows  that  field  enhancement  can  be  great  at  the  points  of  high  curva¬ 
ture  (the  so-called  lightning  rod  effect).3  However,  the  most  effective  SERS  systems  are 
collections  of  interacting  particles.  The  reason  for  this  is  understood  to  be  even  greater 
field  enhancement  that  results  from  electromagnetic  interactions  among  the  constituent 
small  metal  particles.  Models  have  been  reported  which  account  for  these  interactions 
in  terms  of  the  modifications  of  the  Lorentz  local  field  operating  on  a  particle  due  to  the 
presence  of  neighboring  particles4.  However,  the  mean  concentration  of  metal  needed 
in  these  treatments  to  achieve  the  red  shifts  observed  experimentally  is  much  too  high. 
Additionally,  Aravind  et  al.  have  calculated5  the  magnitude  of  the  electromagnetic 
field  at  various  locations  in  the  vicinity  of  two  interacting  particles  invoking  very  high 
multipoles.  Also,  Liver  and  others  have  reported6  the  results  of  a  calculation  among 
several  metal  particles,  but  restricting  the  interaction  to  dipole-dipole.  These  calcula¬ 
tions  demonstrate  that,  m  addition  to  a  further  increase  in  the  local  field  magnitude, 
the  interaction  among  particles  also  results  in  a  red  shift  of  the  maximum  enhancement 
wavelength  from  the  surface  plasmon  resonance  positioned  in  the  near  uv.  However, 
none  of  the  above  mentioned  theories  can  satisfactorily  explain  the  form  of  the  spectral 
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contour  of  SERS,  which  is  very  wide  with  the  maximum  enhancement  in  the  red  region. 

Although  the  precise  structure  of  roughened  metal  surfaces  is  not  known  with 
certainty,  the  structure  of  colloid  aggregates  has  been  investigated  by  various  methods 
and  shown  to  be  fractal  (see  e.g.  the  electron  microscopy  study  of  Ref.  7).  The  Hausdorff 
dimension  D  of  such  aggregates  is  consistent  with  that  characteristic  of  cluster-cluster 
aggregation7-9 .  For  the  sake  of  brevity,  we  will  refer  to  fractal  clusters  simply  as  fractals. 
This  should  not  cause  confusion  because  no  other  physical  realizations  of  fractals  will 
be  discussed. 

A  tentative  theory  of  enhanced  Raman  scattering  by  fractals  has  already  been 
considered.10  Because  it  is  based  on  a  binary  approximation,  this  treatment  is  only 
applicable  to  diluted  fractals,  making  it  unable  to  account  for  the  form  of  either  the 
absorption  spectrum  or  the  Raman  excitation  contour.  A  scale  invariant  theory  has 
been  developed  in  Ref.  11  to  describe  the  linear  optical  polarizability  of  fractals.  It 
demonstrates  the  scaling  properties  of  the  absorption  and  determines  the  dispersion 
relation  for  the  dipolar  excitation  modes  of  the  fractal  (see  also  Ref.  12).  The  important 
qualitative  feature  of  scaling  are  the  dominance  of  the  strong  fluctuations  of  local  fields 
at  all  scales  and  a  dramatic  restriction  in  the  spatial  extent  of  the  eigenmodes  to  small 
regions  within  the  cluster.  Such  localization  for  the  fractal  vibrations  (“fractons”)  has 
been  predicted13  by  Alexander  and  Orbach,  and  experimentally  observed14  recently  by 
Sapoval  and  others  for  acoustic  excitation  of  fractal  drums. 

A  scale-invariant  theory  developed  by  S.  Alexander15  describes  Raman  scatter¬ 
ing  from  vibration  modes  of  fractals.  However,  this  theory  does  not  include  dipolar 
excitations  (plasmons)  and,  therefore,  can  not  describe  the  phenomenon  of  SERS.  In 
the  present  paper,  we  develop  a  scale-invariant  theory  of  SERS  from  fractals  which 
is  insensitive  to  details  of  the  cluster  structure  providing  universal  description  of  this 
phenomenon  in  terms  of  a  few  critical  exponents.  The  present  theory  fully  takes  into 
account  the  strong  spatial  fluctuations  of  local  fields  characteristic  of  fractals  and  their 
excitations.  It  is  capable  of  explaining  both  the  magnitude  and  spectral  contour  of  the 
enhancement. 
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In  Sec.  2  we  formulate  the  model  and  basic  equations  governing  linear  response 
of  the  fractal.  In  Sec.  3  general  expressions  for  the  enhancement  factor  GRS  are  obtained, 
and  self-similarity  is  invoked  to  derive  the  scaling  expressions  for  GRS .  Numerical 
simulation  with  the  use  of  the  Monte  Carlo  method  in  Sec.  4  confirms  the  predictions 
of  the  scaling  theory.  The  results  are  discussed  and  shown  to  be  in  a  good  agreement 
with  experimental  data  in  Sec.  5. 

2.  FORMULATION  OF  THE  MODEL  AND  LINEAR  RESPONSE  OF  A 
CLUSTER 

Let  us  very  briefly  recapitulate  the  formulation  of  the  model  and  the  results  for 
the  linear  optical  response  of  fractals  (Ref.  11),  which  are  necessary  for  the  present  paper. 
A  fractal  is  considered  as  a  system  of  N  polarizable  particles  (monomers)  with  dipole- 
dipole  interaction  between  them  at  optical  frequencies.  The  monomers  are  positioned 
at  points  r ,,  i  =  1, . . . ,  N.  The  number  iV  of  monomers  in  a  fractal  scales  as 

N  ~  (RC/R„)D  ,  (1) 

where  D  is  an  index  called  the  Hausdorff  dimension,  Rc  is  a  characteristic  total  size  of 
the  fractal,  and  Rq  is  a  constant  equal  to  a  typical  separation  between  monomers.  A 
fractal  is  called  nontrivial  if  D  <  3.  The  length  Rq  plays  the  role  of  a  minimum  scale 
of  the  fractal,  and  the  dependence  on  Rq  is  of  principal  importance  to  the  theory  (see 
Ref.  11  and  below). 

Since  Raman  scattering  is  accompanied  by  a  Stokes  frequency  shift,  the  linear 
response  of  the  cluster  at  the  frequency  uj  of  the  exciting  radiation  can  be  found  inde¬ 
pendently  from  its  Raman  response  at  the  shifted  frequency  u>s.  Experimentally,  most 
fractal  clusters  are  larger  than  the  exciting  wavelength  A.  However,  the  localization  (co¬ 
herence)  length  of  fractal  excitations  can  be  much  smaller  them  A.  Also,  under  certain 
conditions,  the  interaction  of  monomers  at  distances  greater  than  A  can  be  neglected. 
We  will  find  qualitative  conditions  for  the  validity  of  the  two  assumptions  at  the  end  of 
the  next  section  [see  Eqs.  (39)  and  (40)]. 
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Accepting  these  assumptions,  one  can  neglect  the  spatial  variation  of  the  electric 
field  E*0)  of  the  exciting  optical  wave  and  obtain  the  well-known  system  of  equations 
describing  the  transitional  dipole  moment  d,,  induced  on  the  tth  monomer  (oscillating 
with  frequency  u), 

Zdia  =  E^-f^(ia\W\j^)dia  ,  (2) 

j=  1 

where  z  =  1, . . . ,  N.  Here  Z  =  Xq1  ,  and  xo  is  the  linear  polarizability  of  a  monomer  at 
the  light  frequency  u>,  and  W  is  the  dipole  interaction  tensor, 

(; ia\W\j(3 )  =  (  r‘>3>  (3) 

1 0,  i=j- 

The  Greek  indices  label  tensor  components,  summation  over  repeated  Greek  indices 
is  implied,  and  rtJ  =  r,  —  ry  and  =  r ,y/r,>.  Of  course,  Eq.  (2)  is  quantitatively 
applicable  to  small  clusters  ( Rc  <))  without  invoking  additional  assumptions  discussed 
above. 

We  shall  reformulate  the  system  (2)  as  an  equation  in  3iV-dimensional  linear 
space,  introducing  the  vector  |d)  with  components  (ia|d)  =  dtQ,  and  similarly  for  other 
vectors.  In  this  way,  we  obtain 

Z\d)  =  \Ew)-W\d)  ,  (4) 

where  W  is  the  operator  determined  by  its  matrix  elements  (3). 

Let  us  introduce  the  eigenvectors  |n)  of  the  W  operator  and  the  corresponding 
eigenvalues  wn:  (n|W|m)  =  wn6nm,  where  n, m  =  1,2, ...,3 N.  Note  that  the  eigen¬ 
values  wn  and  components  of  the  eigenvectors  (ior|n)  axe  all  real  due  to  the  symmetry 
of  W.  The  solution  of  Eq.  (4)  in  the  |n)  basis  has  the  form 

(n|d)  =  (n|E(0))An  ,  A„  =  (Z  +  wn)~' .  (5) 

From  Eq.  (5),  the  expression  of  Ref.  11  for  the  polarizability  tensor  x«  of  the  ith  monomer 
in  the  cluster  can  be  obtained  as16 

Xi.ap  =  An(t'a|n)0'/?|n)  .  (6) 

nj 
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The  local  field  Ex  acting  on  the  ith  monomer  can  easily  be  expressed  as 


Eia  =  Xo  lXi,a$&0  •  (7) 

After  averaging  over  the  orientation  of  a  cluster  as  a  whole,  the  polarizability  tensor  is 


reduced  to  a  scalar 


X  =  3  {Xi,oa)  , 


where  ('•')  indicates  averaging  over  the  ensemble  of  clusters.  The  total  absorption  cross 
section  <ra  of  the  cluster  is  proportional  to  Imx,  namely  <ra  =  ^irkNlmx,  where  k  is 
the  light  wavevector.  Below  we  will  refer  to  Imx  simply  as  the  absorption. 

The  exciting  light  frequency  u>  enters  Eq.  (4)  and  subsequent  expressions  im¬ 
plicitly,  via  the  complex  variable  Z.  Also,  ail  material  and  geometrical  properties  of 
monomers  affect  the  problem  only  through  Z.  Let  us  isolate  the  real  and  imaginary 
parts  of  Z  in  the  form  X  =  — ReZ,  8  =  —  ImZ.  The  quantity  X  plays  the  role  of  a 
spectral  variable  in  place  of  and  the  parameter  8  >  0  describes  the  dissipation  in  a 
monomer.  The  dependence  of  both  X  and  8  on  u;  for  real  systems  is  discussed  in  Sec.  5. 

Another  quantity  of  fundamental  importance  for  the  theory  is  the  density  v(X) 
of  fractal  eigenmodes  per  monomer  defined  as 

u(X)  =  N~l(s28(X-wn)J  •  (9) 

The  absorption  and  eigenmode  density  satisfy  the  exact  sum  rules11 


lmX(X)dX  =tt  ,  f  u(X)dX  =  3 

J — oo  •/ — oo 


The  principal  advantage  of  using  X  is  shown  in  Ref.  11  to  be  due  to  the  fact 
that  optical  responses  scale  with  X  but  not  with  u.  In  particular,  Imx(^Y)  and  X ) 
have  identical  scaling, 


lmx(X)  ~  WX)  ~  W)'-  '  .  (11) 

Here  d0  is  a  critical  index  introduced  in  Ref.  11,  which  governs  the  properties  of  optical 
responses  of  fractals  and  is  called  the  optical  spectral  dimension.  It  plays  a  similar  role 
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for  polar  excitations  of  fractals  as  Alexander  and  Orbach’s  index13,  fracton  dimension 
d ,  does  for  vibrational  excitation.  However,  in  contrast  to  d ,  the  physical  range  for  da 
is  0  <  d„  <  1.  The  condition  of  scaling  has  the  form11 

6,  <  |jf|  <  R ~3  .  (12) 

For  X  in  this  region,  the  light  excites  collective  modes  of  the  fractal,  i.e.  modes  delo¬ 
calized  over  many  monomers.  However,  these  modes  axe  still  well  localized  within  the 
whole  fractal.  In  other  words,  the  coherence  length  Lx  of  the  excitations11, 

Lx  ~  i*o  (*o  1*!)"^  ,  (13) 

conforms  to  the  inequality  Rq  <C  Lx  Rc •  Due  to  this,  the  optical  responses  of  the 
fractal  per  monomer  do  not  depend  either  on  the  fractal  fine  structure  or  on  the  total 
number  of  monomers  N. 

Since  there  is  no  dependence  on  the  fractal  fine  structure  in  the  range  given 
by  Eq.  (12),  we  can  always  unify  some  number  of  nearest  monomers  into  a  composite 
particle  and  call  it  a  new  (renormalized)  monomer.  Such  a  renormalization  transforma¬ 
tion  changes  only  the  fractal  structure  on  a  minimum  scale;  in  particular,  the  length 
Ro  is  changed.  The  total  absorption  of  the  fractal  should  not  change  under  the  renor¬ 
malization  transformation,  which  means  that  cra  oc  R%.  From  this  and  Eq.  (11),  a 
transformation  law11  follows, 

\X\  oc  .  (14) 

3.  ENHANCED  RAMAN  SCATTERING  BY  CLUSTERS 

We  assume  that  each  monomer  of  the  cluster  apart  from  the  linear  polarizability 
Xo  possesses  also  Raman  polarizability  «.  The  exciting  field  E  applied  to  an  isolated 
monomer,  therefore,  induces  a  dynamic  dipole  moment  dJ  oscillating  with  the  Stokes- 
shifted  frequency  us.  In  order  to  avoid  unnecessary  complications,  k  will  be  assumed  to 
be  scalar.  Accordingly,  we  have  da  =  kE.  The  Raman  polarizability  may  be  either  the 
polarizability  of  a  monomer  itself  or  that  of  adsorbed  molecules  bound  to  the  monomer. 
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For  spontaneous  Raman  scattering,  which  is  an  incoherent  optical  process,  the 
Raman  polarizabilities  /e,  of  different  monomers  (t  =  1,2,... , N)  contain  un correlated 
random  phases.  This  can  be  expressed  as 

(«?«,•)=  W2  Sij  ■  (15) 

This  feature  constitutes  the  principal  difference  between  k  and  the  linear  polarizability 
Xoi  ensuring  that  there  will  be  no  interference  between  the  Stokes  waves  generated  by 
different  monomers. 

The  field  acting  upon  an  zth  monomer  in  a  cluster  is  the  local  field  E,  rather 
than  the  external  field  E*°\  Likewise,  it  is  the  dipole  interactions  of  the  monomers  at 
the  Stokes-shifted  frequency  u>a  that  should  be  considered.  The  components  of  d'  obey, 
therefore,  the  following  system  of  equations 

-  xi  ,  (16) 

where  \o  is  the  linear  polarizability  of  an  isolated  monomer  at  the  Stokes-shifted  fre¬ 
quency  u,.  Equations  (16)  can  be  rewritten  as  a  vector  equation  in  31V-dimensional 
space 

\<f)  =  k\E)  -  xlW\d‘)  ,  (17) 

where  k  is  an  operator  defined  by  its  matrix  elements 

(t'o|/£|j^)  =  KiSijSo0  .  (18) 


The  formal  solution  to  Eq.  (17)  has  the  form 

I"')  =  iriw  • 

where  Z$  =  (xo)_1-  Combining  Eqs.  (4)  and  (19),  one  obtains 


IO  = 


z,  +  w  z  +  w 


|E<°>) 


Using  the  completeness  of  vector  sets  |n)  and  |i),  we  rewrite  (20)  in  the  form 
d'ia  =  Z*Z  KiK&n'  (iot|n)  (j/?|n)  (;7?|n')  (j'0'\ n)  E(0V  . 


(19) 

(20) 

(21) 
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Using  Eq.  (6)  to  clear  the  matrix  elements  from  (21),  we  obtain  the  following  expressions 
for  the  total  Stokes  dipole  moment  Ds  of  the  cluster 


=  =  ,  (22) 

«  J 

where  =  \j(Xa),  as  given  by  Eq.  (6)  with  substitution  of  Xs  =  — R eZa  for  X ,  is  the 
polarizability  of  the  ;  th  monomer  in  the  cluster  at  frequency  ujs. 

The  enhancement  factor  of  SERS  is  given  by 


GRS  = 


(|D.P) 

N  |k|2  |E(°)f 


This  factor  expresses  increase  in  the  Raman  intensity  of  light  scattered  by  a  fractal 
cluster  consisting  of  N  monomers  with  respect  to  the  integral  Raman  intensity  of  light 
scattered  by  N  separate  (non-aggregated)  monomers. 

Expressions  (16)-(23)  are  valid  for  clusters  of  any  geometry,  fractal  or  not.  To 
find  a  closed-form  expression  for  GRS ,  we  must  invoke  the  concept  of  scaling  as  well  as 
some  approximations.  First,  we  consider  the  case  of  the  Stokes  shift  so  large  that  the 
Raman-scattered  light  is  well  out  of  the  absorption  band  of  the  fractal.  This  requires, 
in  particular,  that  |JV  —  Xa\  1.  In  this  case,  the  polarizability  (6)  at  frequency 
is  expressed  as  x*  Qg  ~  Z~x8ag.  Substituting  this  into  Eq.  (22),  after  averaging  over  all 
orientations,  we  obtain  from  Eq.  (23) 


The  exact  relation  Yi  Xi,agXi,ay  —  <5-1  has  been  proved  in  Ref.  11,  based 

on  the  completeness  of  the  eigenvector  set  (iajn).  Substituting  this  into  Eq.  (24)  and 
taking  account  of  Eq.  ^7),  we  finally  obtain 

Grs  =<5(1  +  X2/S2)lmx  •  (25) 

This  expression  for  GRS  coincides  exactly  with  the  expression  found  previously  [see 
Eq.  (20)  in  Ref.  11]  for  the  enhancement  factor  G  of  the  local  field  fluctuations.  This  is 


9 


expected,  because  in  the  limiting  case  considered,  the  Raman-scattered  light  does  not 
interact  with  the  fractal,  and  the  Raman  scattering  intensity  is  simply  proportional  to 
the  mean  square  of  the  local  fields.  The  exact  expression  (25)  reduces  in  the  scaling 
region  to  the  form 

G*5~X2*-IImx~(?(i$|X|)i*+1  ,  (26) 

where  Q  =  (RqS)  1  is  a  large  dimensionless  parameter  [cf.  Eq.  (12)]  which  plays  the 
role  of  the  resonance  quality  factor  for  the  monomer. 


If  the  Stokes  shift  is  not  very  large,  the  general  expression  for  GRS  is  needed, 
which  follows  from  Eqs.  (6),  (22)  and  (23)  after  averaging  over  the  orientations, 


(27) 


where  the  kernel  K  is  defined  as 


Knn.  mm'  —  anaan'  0amaam‘  0  Y  CHn)  (nK)  0^lm)  O'^K)  (28) 

j 

with  ana  =  («*|rc).  The  completeness  and  orthogonality  of  the  eigenvectors  |n) 

ensures  that  the  kernel  K  satisfies  the  sum  rule 

Y  K™‘mm>  =  3  N  .  (29) 

nn'  mm1 

Expression  (27)  can  not  be  evaluated  analytically  in  a  general  form.  Instead, 
we  will  consider  the  limiting  case  of  an  extremely  small  Stokes  shift,  |X  —  X$\  6, 

which  is  experimentally  satisfied  for  most  Raman  bands.  Then  the  Stokes  shift  can  be 
neglected  in  Eq.  (27),  and  one  cam  set  Z  =  Z„  A  =  A*.  In  the  scaling  region  (12),  we 
can  put  |Z|  =  |X|  to  a  high  degree  of  accuracy.  Taking  this  into  account,  we  obtain 
from  Eq.  (27) 

GRS=-^X“g(X)  .  g(X)s  Y,  .  (30) 

nn1  mm' 

To  further  evaluate  GRS ,  we  will  invoke  scale  invariance  and  employ  the  sum 
rule  approach.  An  exact  sum  rule  for  the  function  gr(X)  (30)  has  the  form 

g(X)  dX  =  -27rIm  Y  1 mm '  RnmRn'm*  (Unm'  "b  Rn'm)  ,  (31) 

nn'mm' 
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where 


Rnm  =  {wn  -  wm  -  2 iS)  1  .  (32) 

The  sum  rules  (10)  axe  satisfied,  on  the  order  of  magnitude,  if  we  use  Imx(J0 
and  u(X)  in  the  settling  approximation  (11)  and  expand  the  integrals  over  the  scaling 
region,  |X|  <  Rq3  .  This  means  that  the  number  of  eigenmodes  contained  in  the  scaling 
region  (12)  is  on  the  order  of  the  total  number  3 N  of  eigenmodes.  Because  the  integral  in 
Eq.  (31)  is  rapidly  converging,  it  is  plausible  that  the  integration  over  the  scaling  region 
gives  the  major  contribution  to  it.  In  that  case,  f  g(X)dX  would  be  scale-invariant. 
This  is  what  we  will  assume.  The  scale  invariance  assumption  is  strongly  supported 
by  numerical  results  (see  below).  The  integral  f  g(X)dX  must  be  proportional  to  N, 
in  order  to  yield  GRS  independent  of  N  [see  Eq.  (30)].  Therefore,  its  scale  invariance 
means  that  N~1  f  g(X)dX  depends  on  neither  Rq  nor  Rc  and  can  be  expressed  in 
terms  of  8  only.  Taking  the  dimensionality  into  account,  we  obtain  the  simple  estimate 
/  g{X)dX  ~  N6~3.  Comparing  this  relation  with  Eq.  (31),  we  arrive  at  the  conclusion 
that  in  Eq.  (32)  |u;m  —  w„\  <C  8  in  the  essential  region  of  integration  in  Eq.  (31).  In 
this  case,  neglecting  the  difference  of  eigenvalues  in  Eq.  (32)  and  taking  Eq.  (29)  into 
account,  we  obtain  from  Eq.  (31)  the  sum  rule 

.  (33) 

Considering  the  dimensionality  arguments,  the  scaling  expression  for  g(X)  fol¬ 
lows  from  Eq.  (33) 

g(X)  ~  Nit'S-3  (flj  ffl/  ,  (34) 

where  9  is  am  index.  The  scale  invariance  of  g(X)  implies  that  under  the  renormalization 
transformation  g(X )  a  f?Q.  From  this  requirement  and  taking  into  account  transforma¬ 
tion  law  (14)  and  Eq.  (1),  we  immediately  obtain  8  =  d0  —  1.  This  gives  us  the  scaling 
form 

GRS  ~  Q3  (R3  |Xj)‘io+3  .  (35) 

The  coefficient  of  the  asymptotic  dependence  in  Eq.  (35)  can  be  estimated  using 
sum  rules  (10)  and  (33).  Taking  account  of  the  scaling  relation  (11),  we  can  rewrite 
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Eq.  (34)  as  g(X)  ~  3Imx(AT).  As  comparison  with  Eq.  (10)  shows,  sum  rule 

(33)  is  satisfied  if  the  coefficient  in  this  relation  is  chosen  to  equal  yielding 


i<?3(U3.Y)4Im*(.r)  . 


Because  relations  (35)  and  (36)  are  the  final  results  and  of  principal  importance, 
let  us  indicate  another  way  of  obtaining  them.  The  product  of  the  A  factors  in  Eq.  (27) 
can  be  rewritten  identically  as 

AnA„<AmAm,  =  RnmRnimi  [A„A„<  +  A^A^,,  +  Rnm‘  (A„  —  A^,,)  +  R*m „<  (A^,  —  An*)]  . 

(37) 

To  yield  the  scaling  form  of  GRS ,  the  difference  of  the  eigenvalues  in  Eq.  (37)  should 
be  negligible  with  respect  to  8.  Assuming  scaling  and  neglecting  this  difference,  as  was 
done  above,  and  using  the  orthonormality  relation  of  the  eigenvectors,  we  obtain  from 
Eq.  (2?) 

G*s*6ihtx,£aLlmA"  ■  (38) 

n 

Taking  Eqs.  (6)  and  (8)  into  account,  we  see  that  Eq.  (38)  is  equivalent  to  Eq.  (36). 

Let  us  emphasize  that  in  the  case  of  a  large  Stokes  shift,  expression  (26)  for 
GRS  is  asymptotically  exact,  i.e.  both  the  scaling  index  da  +  1  and  the  proportionality 
coefficient  axe  correct.  In  the  case  of  a  small  Stokes  shift,  we  believe  that  the  index 
d0  +  3  in  Eq.  (35)  is  also  asymptotically  exact  because  it  is  derived  from  the  general 
requirements  of  scale  invariance.  Unlike  it,  the  numerical  coefficient  \  in  Eq.  (36)  is 
only  an  estimate,  its  accuracy  depending  on  what  fractions  of  the  sum  rules  (10)  and 
(33)  are  saturated  in  the  scaling  region  (12). 

The  scale-invariant  theory  presented  above  is  developed  based  on  Eq.  (2),  where 
the  spatial  variation  of  the  exciting  field  E^0^  is  neglected.  A  necessary  condition  for 
the  validity  of  such  neglect  and  also  for  the  existence  of  scaling  is  Lx  A,  which,  upon 
taking  account  of  Eq.  (13),  acquires  the  form 
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This  condition  does  not  contradict  Eq.  (12),  because  it  is  assumed  that  Rq  «C  A,  as  is 
the  case  experimentally  (note  that  D  <  3  and  d0  <  1). 


To  have  the  indices  of  GRS  coincide  with  the  values  d0  +  1  and  d0  +  3  [Eqs.  (26) 
and  (35)],  one  more  condition  must  be  satisfied,  namely  that  the  external  field  at  each 
monomer  is  the  same,  which  is  implied  in  Eq.  (2).  This  assumption  is  rigorously  correct 
only  for  small  clusters  ( Rc  A).  The  present  theory  can  be  applicable  to  the  clusters 
with  Rc  »  A  only  if  the  fax-field  zone  [i.e.  monomers  positioned  at  distances  rtJ  A 
from  a  given  ith  one]  contributes  negligibly  to  the  local  field  E,.  The  condition  for  this 
has  been  found  in  Ref.  17  to  be 


X\  > 


(-Ro/A)3-D , 


for  D  <  2  ; 
for  D  >  2  . 


(40) 


Note  that  for  the  most  common  case  of  D  <  2,  the  first  of  conditions  (40)  supersedes 
(39).  If  this  is  uist,  then  the  far  field  zone  (which  might  otherwise  bring  about  inter¬ 
ference  phenomena)  can  be  neglected.  In  the  intermediate  zone  (r,y  ~  A)  the  external 
field  E(0)  is,  strictly  speaking,  not  constant,  but  does  not  change  phase  often  enough  to 
cause  destructive  interference.  Therefore,  in  the  zeroth-order  approximation,  we  may 
neglect  the  variation  of  E*0)  and  use  Eq.  (2). 


Let  us  briefly  recapitulate  and  discuss  the  analytical  results  obtained  above. 
The  enhancement  factor  of  Raman  scattering  GRS  for  each  of  the  cases  of  large  and 
small  Stoke  shifts  [see  Eqs.  (26)  and  (35),  (36),  respectively]  has  the  form  of  a  power 
law  in  the  variable  AC,  with  the  corresponding  index  determined  by  the  optical  spectral 
dimension  d„,  provided  X  is  within  the  region  given  by  Eq.  (12).  For  the  case  of  large 
shifts,  this  result  is  rigorously  shown  to  be  asymptotically  exact,  with  an  index  equal 
1  +  d0.  For  small  Stokes  shifts,  scale  invariance  of  GRS  is  assumed  [this  assumption 
is  strongly  supported  by  the  results  of  numerical  simulation  (see  below)],  and  then  the 
index  3  4-  d0  is  obtained.  Both  the  indices  mentioned  above  are  positive,  implying  that 
the  enhancement  increases  with  X.  Physically,  this  property  follows  from  the  decrease 
of  the  eigenmode  coherence  length  Lx  [see  Eq.  (13)],  which  brings  about  enhancement 
of  the  local  field  fluctuations  (from  one  monomer  to  another)  in  the  fractal. 
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The  Raman  scattering  enhancement  factor  GRS  is  found  to  be  large,  propor¬ 
tional  to  large  factors,  namely  Q  for  the  case  of  small  and  Q3  for  large  Stokes  shifts. 
Note  that  Q  is  required  to  be  a  large  quantity  as  a  necessary  condition  for  the  very 
existence  of  the  scaling  region  [see  Eq.  (12)].  Thus,  strong  enhancement  of  Raman  scat¬ 
tering  is  characteristic  of  scaling  and,  consequently,  of  fractals  consisting  of  monomers 
with  high-quality  optical  resonances  (i.e.  with  Q  1),  such  as  silver  and  gold  fractal 
clusters. 

The  Raman  scattering  intensity  is  proportional  to  the  intensity  (mean-squared 
magnitude)  of  the  local  fields  induced  on  different  monomers.  Strong  fluctuations  of  the 
local  fields  in  fractals  bring  about  an  increased  local  field  intensity  and,  consequently, 
enhanced  Raman  scattering.  In  a  pure  form  this  mechanism  works  for  large  Stokes 
shifts,  in  which  case  the  scattered  light  does  not  interact  with  the  cluster.  Therefore, 
the  corresponding  expression  (25)  GRS  simply  coincides  with  the  coefficient  G  of  the 
enhancement  of  local  fields  found  in  Ref.  11. 

Apart  from  the  general  mechanism  discussed  above,  for  the  case  of  a  small  Stokes 
shift  there  exists  an  additional  enhancement.  The  scattered  field  is  not  radiated  freely, 
but  rather  interacts  with  the  cluster.  More  exactly,  the  dipole  moment  dJ,  induced  at 
the  shifted  frequency  u,,  induces,  in  turn,  secondary  dipole  moments  at  u3,  and  so  on 
to  self-consistency.  The  Raman  radiation  is  emitted  by  the  integral  dipole  moment  D' 
(22),  resulting  in  an  additional  enhancement  as  given  by  Eq.  (35).  We  point  out  that 
for  the  case  of  a  small  Stokes  shift,  GRS  (35)  is  certainly  not  the  square  of  GRS  given 
by  Eq.  (26)  for  large  Stokes  shifts,  as  is  often  proposed  in  the  SERS  literature. 

4.  NUMERICAL  MODELING 

We  have  examined  two  types  of  fractals,  the  random  walk  fractals  [fractal  di¬ 
mension  D  =  2  and  optical  spectral  dimension  d0  a;  0.4  (Refs.  11  and  17)]  and  cluster- 
cluster  aggregates8  [D  «  1.7  (Refs.  7  and  9)  and  d0  as  0.3  (Ref.  17)].  The  clusters  were 
generated  using  the  Monte  Carlo  method  and  well-known  procedures.  Then,  following 
Ref.  11,  the  clusters  have  been  subjected  to  dilution  (decimation),  which  consists  of 
the  following.  The  ith  (i  =  1,2,...,  N)  monomer  is  randomly  retained  in  the  cluster 
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with  some  probability  d,  or  removed  with  probability  1  —  3.  This  procedure  simplifies 
the  fractal  structure  at  small  scales  and  reduces  the  total  number  N  of  monomers  in 
the  cluster  on  average  by  a  factor  simplifying  greatly  the  numerical  calculations.  At 
the  same  time,  the  resulting  (diluted)  fractal  is  characterized  by  the  same  HausdorfF 
dimension  as  the  original  one.  In  most  cases,  32-fold  decimation  (/?  as  0.03)  has  been 
performed.  Some  simulations  have  been  done  with  3  as  small  as  10~3  for  comparison. 
The  results  of  the  computations  clearly  show  that  the  optical  properties  of  fractals  in  the 
scaling  region  [i.e.  for  X  satisfying  Eq.  (12)]  do  not  depend  on  the  dilution,  as  expected. 
Finally,  averaging  over  a  large  ensemble  of  fractals  (up  to  103)  has  been  performed. 

The  enhancement  factor  GRS  of  the  cluster-enhanced  Raman  scattering 
(CERC)  is  calculated  accordingly  to  the  exact  expression  (23)  with  the  use  of  Eqs.  (6) 
and  (22).  Most  interesting  from  the  experimental  point  of  view  is  the  case  of  small 
Stokes  shifts,  where  GRS  is  greatest.  Also,  because  expression  (26)  for  GRS  for  the  case 
of  large  Stokes  shift  is  asymptotically  exact,  its  properties  are  clem.  It  is  for  this  reason 
that  we  concentrate  below  on  the  numerical  modeling  for  the  case  of  small  Stokes  shifts. 

The  theory  [Eqs.  (35)  or  (36)]  predicts:  (i)  scaling  of  GRS  as  a  function  of  the 
variable  X  with  index  d0+ 3;  ( ii )  proportionality  of  GRS  to  Q3;  (Hi)  independence  of  GRS 
from  jV;  and  (iv)  equality  of  the  scaling  indices  of  GRS  for  X  >0  and  X  <  0.  The  two 
last  features  are  not  self-evident.  In  fact,  there  are  four  summations  over  n  =  1,2, ...,  ZN 
in  the  exact  expression  (30),  and  it  is  not  clear  a  priori  that  these  summations  mutually 
cancel  to  result  in  the  independence  of  GRS  from  N.  The  equality  of  indices  for  X  >  0 
and  X  <  0  is  not  evident  because  Eqs.  (2)  and  (16),  along  with  all  the  solutions  which 
follow  from  them,  axe  not  even  with  respect  to  the  sign  change  X  <-►  —X.  The  numerical 
computations  strongly  support  these  analytical  results. 

The  normalized  enhancement  factor  GRSQ~3  for  dusters  generated  assuming 
the  cluster-cluster  aggregation  is  shown  in  Fig.  1  as  a  function  of  the  dimensionless 
spectral  variable  | JY  |  for  both  positive  and  negative  X  (note  the  double-logarithmic 
scale).  The  calculated  points  lie  along  straight  lines,  indicating  scaling  of  GRS  over 
seven  decades  of  its  magnitude.  The  scaling  indices  found  from  these  data,  3.30  for 
X  >  0  and  3.27  for  X  <  0,  are  almost  equal  and  close  to  the  predicted  value  of  da  +  3 
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with  the  optical  spectral  dimension  found  in  Ref.  17  of  d0  ~  0.3.  The  lines  for  X  >  0 
and  X  <  0  are  parallel  but  not  coincident,  which  shows  that  the  coefficients,  unlike  the 
indices,  are  unequal  (they  differ  nearly  by  a  factor  of  2).  The  last  property  results  from 
the  lack  of  symmetry  with  respect  to  the  sign  change  of  X,  emphasizing  the  nontriviality 
of  the  equality  of  indices  for  X  >  0  and  X  <  0. 

Further  insight  into  the  requirements  for  scaling  and  the  test  of  the  proportion' 
ality  GRS  oc  Q3  is  provided  by  Fig.  2,  where  GRSQ~3  is  plotted  as  a  function  of  R%X 
for  two  different  values  of  Q.  We  see  that  for  Q  =  200  scaling  takes  place  in  the  whole 
region  presented  in  Fig.  2,  while  for  Q  =  50  the  scaling  exists  only  when  R%X  >  0.1, 
where  the  two  sets  of  data  virtually  coincide.  This  shows  the  importance  of  the  scaling 
condition  (12),  which  requires  XR%  »  Q~l  The  independence  of  GRSQ~3  from  Q 
in  the  region  of  scaling,  clearly  seen  in  Fig.  2,  corroborates  the  analytical  result  that 
G  oc  Q3  in,  and  only  in,  this  region. 

Comparison  of  the  results  obtained  from  approximate  formula  (36)  with  the 
computations  based  on  the  exact  expression  (23)  is  presented  in  Fig.  3.  We  can  see  that 
the  predictions  of  Eqs.  (23)  and  (36)  do  not  differ  by  more  than  a  factor  of  two  over 
seven  decades,  which  makes  Eq.  (36)  a  useful  approximation.  The  indices  of  the  two 
sets  of  data  are  indeed  very  close,  and  a  very  good  quantitative  agreement  is  obtained 
if  the  numerical  coefficient  in  Eq.  (36)  is  changed  from  |  to  j.  The  possible  origin  of 
the  difference  between  the  numerical  coefficients  is  discussed  at  the  end  of  the  preceding 
section. 


The  data  presented  in  Fig.  4  demonstrate  scaling  for  random  walk  fractals,  and 
also  provide  a  direct  test  of  independence  of  the  enhancement  coefficient  GRS ,  as  defined 
by  Eq.  (23),  from  cluster  size.  The  latter  property  is  evident  from  Fig.  4,  where  the 
points  corresponding  to  N  =  64  and  =  1 28  coincide  within  statistical  error.  The 
scaling  is  as  good  as  for  the  case  of  Fig.  1,  and  the  index  found.  3.45,  is  reasonably  close 
to  the  predicted  value  of  3  +  d0  with  d0  m  0.4  as  found  in  Refs.  11  and  17. 
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5.  DISCUSSION 


Since  the  analytical  results  obtained  have  been  discussed  at  the  end  of  Sec.  3,  we 
will  mainly  concentrate  our  discussion  on  the  underlying  physics  and  manifestation  of 
the  predicted  critical  behavior  in  observable  optical  phenomena,  including  quantitative 
comparison  with  experimental  data. 

FYactals  are  self-similar  objects,  and  the  general  property  of  such  objects  is  such 
that  fluctuations  in  space  on  sill  scales  axe  of  the  same  relative  magnitude.  Therefore 
it  is  understandable  that  a  mean-field  approach  is  not  applicable,  and  that  the  scaling 
behavior  of  linear  optical  responses  takes  place,  as  established  in  Ref.  11.  The  linear 
optical  absorption  and  density  of  fractal  eigenmodes  scale  with  the  same  index  d„  —  1, 
where  d„  is  called  the  optical  spectral  dimension.  This  scaling  is  not  present  in  terms  of 
the  frequency,  but  rather  in  terms  of  a  “natural”  spectral  variable  X  =  Rexo"1,  where 
Xo  is  the  optical  polarizability  of  a  single  monomer.  The  variable  X  can  be  expressed 
in  terms  of  light  frequency  u  or  wavelength  A,  but  in  a  complicated,  non-scaling  form 
given  below.  The  region  of  scaling  is  determined  by  Eq.  (12)  with  additional  conditions 
(39)  and  (40)  for  large  ( Rc  A)  clusters. 

The  main  goal  of  this  work  is  to  examine  how  the  scaling  properties  of  the 
fractal,  its  eigenmodes  and  linear  responses  manifest  themselves  in  the  Raman  scattering 
of  light.  The  major  result  of  this  paper  is  the  prediction  of  the  scaling  behavior  of  the 
enhancement  factor  GRS  of  the  Raman  scattering.  Before  proceeding,  let  us  recall  that, 
by  definition,  the  factor  GRS  shows  how  many  times  the  Raman  intensity  per  monomer 
is  increased  upon  aggregation  of  the  monomers  into  a  cluster.  The  Raman  susceptibility 
of  a  monomer  may  represent  either  its  own,  inherent  Raman  scattering  ability,  or  that 
of  the  adsorbates  bound  to  it. 

Physically,  the  intensity  of  Raman  scattering  is  proportions!  to  the  averaged 
intensity  of  the  local  fields  acting  on  the  monomers  in  the  cluster.  Such  a  mechanism, 
acting  in  a  pure  form  in  the  case  of  very  large  Stokes  frequency  shifts,  predicts  GRS  (25) 
to  be  equal  to  the  enhancement  factor  G  of  the  local  fields,  which  has  been  calculated11 
exactly  earlier.  In  this  case,  the  scaling  of  GRS  in  X  follows  from  the  scaling  of  the 
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optical  absorption  [Im\(X)].  The  magnitude  of  GRS  is  large,  proportional  to  the  large 
parameter  Q  (quality  factor  of  the  optical  resonance  of  monomers),  and  the  scaling 
index  is  d0  +  1. 

In  the  case  of  a  small  Stokes  shift,  which  is  typical  for  most  experiments  (in¬ 
cluding  those  discussed  below),  the  dynamic  dipoles  oscillating  at  the  Stokes  frequency 
interact  with  the  fractal,  resulting  in  additional  enhancement  of  Raman  scattering.  In 
this  case,  GRS  (35)  is  much  larger,  proportional  to  Q3,  and  the  index  is  d„  +  3.  We 
note  that  the  tentative  theory  of  Raman  scattering  enhancement  based  on  the  binary 
approximation  (Ref.  10)  predicts  correctly  the  dependence  of  GRS  on  Q,  but  fails  to 
give  the  correct  scaling  in  X. 

The  results  of  the  numerical  simulation  presented  above  (Sec.  4)  strongly  sup¬ 
port  the  scaling  of  GRS  in  terms  of  X,  the  independence  of  GRS  from  the  number  of 
monomers  in  the  fractal,  and  the  proportionality  of  GRS  to  Q3.  Given  the  large  val¬ 
ues  of  Q  characteristic  of  metal  fractal  clusters,  the  last  feature  qualitatively  explains 
the  strong  enhancement  of  Raman  scattering  observed  with  aggregated  colloidal  metal 
clusters. 

Now  let  us  examine  how  the  scaling  relations  established  and  discussed  above 
translate  into  experimentally  measurable  spectral  profiles  of  SERS,  i.e.  how  GRS  de¬ 
pends  on  the  wavelength  A.  To  do  this,  we  need  to  know  the  variables  X  and  6  as 
functions  of  A.  Let  us  concentrate  on  metallic  fractal  clusters.  Such  a  cluster  consists 
of  aggregated  spherical  colloidal  particles  as  monomers.  The  polarizability  xo  of  such  a 
monomer  is  given  by  the  familiar  expression 

Xo  =  R-m  (e  —  fo)  (c  +  2eo)  1  ,  (41) 


where  Rm  is  the  radius  of  the  spherule,  and  e  and  e0  are  the  dielectric  permittivity  of 
the  colloidal  material  (metal)  and  embedding  medium  (usually  water).  From  Eq.  (41), 
the  spectral  variable  X  and  the  resonance  quality  factor  Q  =  l/R^S  are  determined  to 
be 


X  = 


|2e  +  eo|2  —  9 
4  |e  —  e0|2 


Q  = 


( Rm  \  3  |e  —  60|2 
\Ro)  3e"e0 


(42) 


18 


where  e"  =  Ime. 


From  Eqs.  (35),  (41)  and  (42),  we  obtain  an  explicit  expression  for  G*5  in  terms 
of  the  metal  dielectric  permittivity  e, 


yRS 


(Ro/Rm) 


3  do 


e  4-  |eo 


ho 


J„+ 3 


|f  -  fo|2li'  (3e"«o)3 


(43) 


Colloidal  metal  particles  are  known  to  be  fractal  clusters  with  the  HausdorfF 
dimension  close  to  that  predicted7,9  by  the  cluster-cluster  aggregation  model,  D  ~  1.7. 
For  such  clusters  we  can  expect17  d0  ~  0.3  and  Ro  ~  Rm.  An  electron  micrograph 
of  a  silver  colloidal  cluster18  is  shown  in  Fig.  5.  The  cluster,  which  contains  ~  103 
metallic  monomers,  was  produced  as  outlined  in  Ref.  19.  One  can  clearly  see  in  Fig.  5 
the  main  qualitative  feature  of  a  fractal  cluster,  namely  its  self-similarity  on  all  spatial 
scales  between  the  minimum  scale  (distance  between  monomers)  and  the  maximum 
scale  (cotal  size  of  the  cluster):  each  fragment  after  magnification  resembles  the  whole 
cluster.  Also,  the  tenuous  structure  of  the  cluster  is  evident:  there  are  cavities  of  all 
sizes  up  to  the  total  size  of  the  cluster,  and  the  larger  the  fragment  that  is  considered, 
the  lower  its  density. 

In  the  red  region  of  visible  light,  the  dielectric  constant  e'  =  Ree  is  known  to 
be  negative  and  large  in  magnitude  for  most  metals,  especially  the  noble  metals,  i.e. 
—e'  >  eo,e"-  In  this  case  GRS  becomes  very  large, 


For  noble  metals,  the  resonance  quality  factor  Q  (42)  can  be  as  high  as  100,  yielding 
the  upper  estimate  GRS  ~  106.  Another  role  of  the  fractality  is  in  the  formation  of  the 
broad  spectral  contour  of  SERS  displaying  the  dramatic  increase  from  the  blue  to  red 
region  where  Q  is  large. 

The  model  with  dipole-dipole  interaction  considered  above  does  not  include 
some  factors  present  in  read  systems.  In  colloidal  clusters,  monomers  are  touching  or 
almost  touching  each  other  (cf.  Fig.  5).  As  a  result  of  this,  there  may  exist  Ohmic 
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conduction  between  the  monomers.  However,  at  high  frequencies,  the  capacitive  con¬ 
ductivity  (which  is  taken  into  account  by  the  dipole  coupling)  will  prevail  over  the 
Ohmic  conductivity,  as  argued  in  Ref.  11.  Another  factor  is  the  short-range  interactions 
of  nearest  monomers,  namely,  exchange  and  higher- multipole  couplings  (the  exchange 
coupling  is  also  responsible  for  conduction  between  the  monomers).  We  believe  that  the 
long-range  dipole  interaction  dominates  in  determining  the  collective  excitation  spec¬ 
trum  of  a  fractal,  i.e.  in  the  scaling  region.  However,  the  higher-multipole  interactions 
may  change  the  local  fields  acting  on  a  monomer  or  on  a  monomer- adsorbed  molecule, 
as  has  been  shown5  for  two  monomers.  This  factor  can  renormalize  the  magnitude  of 
enhancement  factor  GRS ,  but,  probably,  has  lesser  effect  on  its  spectrum.  We  hope  to 
return  to  quantitative  consideration  of  these  factors  elsewhere. 

Let  us  compare  the  theory  with  experimental  results  obtained  with  colloidal 
silver.  Figure  6  shows  the  spectral  dependence  of  X  and  S  for  silver  calculated  with 
optical  constants  adopted  Tom  Ref.  20.  One  can  conclude  from  these  data  that  the 
scaling  region,  where  X  >■  £,  is  A  >  400  nm.  Also,  in  this  region  Q  =  1/R%&  is  large. 

Experimental  SERS  enhancement  data  obtained  in  Ref.  19  are  compared  in 
Fig.  7  with  GRS  values  calculated  with  Eq.  (43)  as  a  function  of  A.  Only  the  spectral 
dependence  of  GRS  is  informative  in  this  figure  since  only  relative  values  of  GRS  are 
reported  in  Ref.  19.  The  experimental  data  presented  in  Fig.  7  are  normalized  by 
putting  Grs  =  15000  at  460  nm,  which  is  a  reasonable  value.  Clearly,  the  present  theory 
accounts  successfully  for  the  dramatic  increase  in  Raman  enhancement  accompanying 
aggregation,  and  for  the  observed  increase  of  GRS  towards  the  red  despite  the  fact  that 
the  enhancement  for  a  single  silver  colloidal  particle  is  expected  to  peak  in  the  near  uv. 
We  emphasize  that,  aside  from  a  single  multiplicable  factor,  the  calculated  enhancement 
factor  contains  no  adjustable  parameters  ( d0  has  previously  been  found  from  modeling 
of  linear  responses). 
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FIGURE  CAPTIONS 


Fig.  1.  The  normalized  enhancement  factor  for  cluster-cluster  aggregates  plotted 
against  the  dimensionless  spectral  parameter  ( Rq  |X|)  in  a  double  logarithmic 
scale.  The  straight  lines  represent  the  best  power- law  fits  to  the  data  (the  fit 
parameters  are  shown  in  the  figure).  The  parameters  used  in  the  computations 
are  0  =  ■—,  Q  =  200,  and  N  =  64  (for  notations  see  text).  Each  point  is 
obtained  by  averaging  over  250  individual  clusters. 

Fig.  2.  The  normalized  enhancement  factor  GRSQ~3  as  a  function  of  R^X  for  two 
values  of  the  resonance  quality  factor,  Q  =  200  and  Q  =  50.  The  rest  of 
parameters  are  the  same  as  in  Fig.  1. 

Fig.  3.  The  normalized  enhancement  factor  GRSQ~3,  computed  with  the  exact  formula 
(23)  (circles)  and  approximation  (36)  (diamonds),  plotted  versus  —R%X  (for 
X  <  0)  in  a  double  logarithmic  scale.  The  straight  lines  correspond  to  the 
power-law  fits  with  the  parameters  indicated  in  the  figure.  The  clusters  and 
parameters  used  in  the  computations  are  the  same  as  in  Fig.  1. 

Fig.  4.  The  normalized  enhancement  factor  GRSQ~3  for  random  walk  fractals  plotted 
versus  R%X  (for  X  >  0)  in  a  double  logarithmic  scale.  The  total  number  of 
monomers  in  a  cluster  is  N  =  128  with  averaging  over  50  individual  clusters 
(circles)  and  N  =  64  with  averaging  over  250  clusters  (stars).  The  straight  line 
corresponds  to  the  best  power-law  fit  with  the  parameters  shown  in  the  figure. 

Fig.  5.  Electron  micrograph  of  a  silver  colloidal  cluster,  adopted  from  Ref.  18.  For 
details  of  obtaining  this,  see  Ref.  19. 

Fig.  6.  Spectral  parameter  X  and  dissipation  parameter  6 ,  both  multiplied  by  R%  to 
give  dimensioneless  quantities,  plotted  against  wavelength  A.  Optical  constants 
for  silver  Eire  taken  from  Ref.  20,  and  the  ratio  Ro/Rm  =  0.7.  typical  for  cluster- 
cluster  aggregation,  is  assumed. 

Fig.  7.  Theoretical  and  experimental  enhancement  factors  for  the  silver  colloid  clusters 
as  functions  of  wavelength  A. 
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Fig.  1  •  M.l  .Stockman,  V.M.Shalaev,  M.Moskovits,  R.Botet,  and  T.F. George. 

Enhanced  Raman  scattering  by  fractal  clusters: 

Scale  invariant  theory 
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Fig. 3.  M.l  .Stockman,  V.M.Shalaev,  M.Moskovits,  R.Botet,  and  T.F. George 

Enhanced  Raman  scattering  by  fractal  clusters: 

Scale  invariant  theory 


Fig. 6.  M.l. Stockman,  V.M.Shalaev,  M.Moskovits,  R.Botet,  and  T.F.George. 
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Fig. 7.  M.l. Stockman,  V.M.Shalaev,  M.Moskovits,  R.Botet,  and  T.F. George. 
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